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Abstract — This paper presents design and performance analy- 
sis of a modified reference model MRAC (M-MRAC) architecture 
for a class of multi-input multi-output uncertain nonlinear 
systems in the presence of bounded disturbances. M-MRAC 
incorporates an error feedback in the reference model defini- 
tion, which allows for fast adaptation without generating high 
frequency oscillations in the control signal, which closely follows 
the certainty equivalent control signal. The benefits of the method 
are demonstrated via a simulation example of an aircraft’s wing 
rock motion. 

I. Introduction 

Controlling uncertain nonlinear systems is a challenging 
task, and remains one of the active research topics in the 
systems theory. There are several directions in this field, one 
of which is the adaptive control. The asymptotic behavior of 
adaptive systems has been a well researched topic during the 
last couple of decades. However, the transient behavior of the 
input and output signals is still a challenging problem. Since 
the transient of the adaptive signals can be very oscillatory 
with big excursions [18], there has been a great deal of effort 
to modify the control architecture and the adaptive laws from 
the perspective of improving it. The majority of these efforts 
led to nonadaptive high gain feedback [3], [4], [16], switching 
control law [9], [10] or to a parameter dependent persistent 
excitation condition [1], and addressed only the behavior of 
output signal. 

First contribution to transient analysis of the the adaptive 
control signal can be found in [6], where it is shown that 
the bound on the control signal is proportional to the square 
root of the adaptation rate. This result is conservative, but 
it reflects the general observations about the control signal 
behavior of the MRAC system. In [2], an adaptive control 
architecture, called L t adaptive control, has been introduced, 
which can achieve arbitrarily close tracking of a given refer- 
ence command both in transient and steady state by increasing 
the adaptation gain. 

In [15] we have introduced the concept of a M-MRAC 
architecture for linear systems that can achieve desired level of 
accuracy in tracking both input and output signals of a given 
reference model by a proper selection of design parameters. 
In this paper we apply the M-MRAC approach to a class of 
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multi-input multi-output uncertain nonlinear systems subject 
to bounded disturbances. The proposed algorithm guarantees 
tracking performance both in input and output signals similar 
to L\ adaptive control, but has a simpler structure and is easier 
to implement. Moreover, it requires the selection of only two 
control parameters, for which a design guideline is provided. 
This guideline is based on the results obtained for second order 
linear time variant systems, which is a contribution by itself. 

The rest of the paper is organized as follows. In Section II 
we state the problem, and in Section III present the control 
design. The error signals are defined in Section IV. In Section 
V the controller’s performance is analyzed and the design 
specifics are discussed. A simulation example is presented in 
Section VI and some concluding remarks are given in Section 
VII. 

II. Problem Formulation 

Consider a multi-input multi-output controllable uncertain 
system 

x(t) = Ax(t) + B[u(t) + W f(x(t)) + d(t)\ (1) 

with x(0) = Xq , where x £ R n and u £ R 9 are the state 
and input of the system, / : R™ — > R 9 is a known vector 
of regressor functions, assumed to satisfy the existence and 
uniqueness conditions, W £ W" /q is a matrix of unknown 
constant parameters, and d : R — > R p is a bounded but 
otherwise unknown external disturbance, and A £ R" x ” and 
B £ R nxp are unknown constant matrices satisfying the 
following matching conditions. 

Assumption 2.1: Given a Hurwitz matrix A rn £ R” x " and 
a matrix B m £ M" xp of full column rank, there exists a matrix 
A'i £ R pxn and a sign definite matrix A £ R pxp such that 
the following equations hold 

B = B m A (2) 

A = A m + BKi . 

Remark 2.1: The sign definiteness of A corresponds to the 
conventional condition on the high frequency gain matrix 
of MIMO systems (see for example [11]). Without loss of 
generality we assume that A is positive definite. The rest of 
the conditions for the existence of an adaptive controller are 
given by the equations (2). 

We notice that systems in the form of (1) frequently arise in 
aerospace applications (see for example [8], [17] for diagonal 
A) and in robotics. Obviously, any fully actuated mechanical 
system can be described by equation (1). 



The objective is to design a control signal u(t) such that 
the state of the system tracks the state x°(t) of a reference 
model 

x°(t) = A m x°(t) + B m r(t), as°(0) = x 0 , (3) 

where A m , B rn are chosen according to performance speci- 
fications and satisfy Assumption 2.1, and r(t) is a bounded 
and piecewise continuous external command. To achieve this 
objective we use the M-MRAC architecture introduced in [15], 
where the system (3) is called an ideal reference model. It is 
not a part of the control design and is only used for the analysis 
purposes. 

Taking into account Assumption 2.1 we write 

x(t) = A m x(t ) + B m r(t ) + B m A[u(t) + 9<f>(t) + d(t)} (4) 

where for the convenience we denote K 2 — —A -1 , <f>(t) = 
[f T (x(t)) x T (t) r T (t)] T , and 9 = [W K\ K 2 \. 

III. Control Design 

According to M-MRAC architecture, the design of the 
adaptive control is based on the modified reference model 

— A rn x rn {t) A B m r(t) -f- Ac(f), ar m ( 0 ) — xq , (5) 

where e(t) = x(t) — x m (t) is the error between the system 
and the modified reference model, A > 0 is a feedback gain 
to be specified in the analysis. The adaptive control is given 
by the equation 

u(t) = — 9{t)4>{t ) — d(t) , 

where 9(t) is the estimate of the unknown matrix 9, and d(t) 
is the estimate of a constant vector d , which is the constant 
part (or an average value) of d(t). The ideal version of this 
control signal that exactly cancels the uncertainties is 

u°(t) = —94>(t) — d(t) . (6) 

Obviously, the ideal control u°(t) reduces the system (4) into 
the ideal reference model (3), which always can be specified 
from the performance perspective. However, the ideal control 
signal ( 6 ) cannot be implemented and is only used for the 
analysis purposes. 

The adaptive laws for the estimates 9(t ) and d(t) are defined 
using the projection based adaptive law 

fa) = 7Pr(<9(i), BlPe(t)<t> T (t)) 

d(t) = 7 Pr (d(f), B^-Pe(f)) , (7) 

where 7 > 0 is the adaptation rate, P = P T > 0 is the 
solution of the Lyapunov equation 

A^P + PA m = -Q ( 8 ) 

for some Q = Q T > 0, and Pr(-,-) denotes the projection 
operator [12] defined as Pr (9,y) = [I — G(9)\y, where 
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with the notation S7w[9) = dip }P , and the smooth convex 
functions ip(9) is given by 


,/h tr (9 t 9) - 6 »^ ax 

= ^2 


with 9 ni;tx denoting the norm bound imposed on the parameter 
matrix 9 and eg denoting the convergence tolerance. The 
projection operator has the following properties 

Lemma 3.1: [12] Let 9 0 G Ll 0 = {9 G R" | <p(9) < 0}, 
and let the parameter 9(t) evolve according to the following 
dynamics 


9(t) = Pr (9(t), y), 9(t 0 ) G fl. (10) 

Then 1) 9(t) G fli = {9 £ R” | c p(9 ) < 1} or ||0(f)|| < 
9* for all t > to, where 9* = y/1 + eg 9 max , 2) [ 9(t ) — 
6» o ] T [Pr(0(f), y) - y] < 0 for all t > t 0 . 

IV. Error Dynamics 

Introducing the parameter estimation error as 9(t) = 9{t) — 
9, the dynamics of the tracking error e(t) can be written in 
the form 

e(t) = (A m - AI„)e(f) + B m A[-9(t)4>(t) + d(t) - d(t)] , 

where I n denotes n-dimensional identity matrix. The error 
signal e(t) is used in stability analysis, but for the performance 
analysis we also need the actual tracking error e°(f) = 
x(t) — x°(t), which satisfies the equation 

e°(f) = A m e°(t) + B m A[-9(t)<p(t) + d(t) - d(t)] . 

These two error signals are related via the linear equation 

^[e(t) - e°(f)] = A m [e(t) - e°(f)] - A e(t) . ( 11 ) 

Since A rn is Hurwitz, the C\ norm of the the state transition 
matrix $(f) = e A ™t j s bounded. That is, there exists a positive 
constants k m such that ||4 > (f )||/; 1 < k rn ■ Therefore, it follows 
from the equation ( 11 ) that 

ll e r(^)IUoo < (! + Afc m )||e r (t )||£ 00 , ( 12 ) 

where the subscript r indicate the extended C 0 0 norm on the 
interval 0 < t < r (see [7], p. 200 for details). Moreover, if 
e(t) G £ 00 . then 

||e 0 (f)||£ oo < (1 + Afe ro )||e(f)|| £oo , (13) 

In the following analysis we will also need the control error 
signal that is defined as u(t) = u(t) — u°(t). From M-MRAC 
architecture it follows that 


u(t) — d(t) — d(t) — 9(t,)4>(t) . (14) 

Therefore, the error dynamics (11) can be also represented as 

e(t) = (A m - AI„)e(f) + B m Au(t ) , (15) 

Since the ideal control signal is the best achievable signal, we 
are interested in minimizing the control error u(t), as well as 
the tracking error signals e(t) and e°(f) by selecting proper 
values for the adaptation rate 7 and feedback parameter A. 



V. Analysis of the M-MRAC Performance 
A. Boundedness 

Theorem 5.1 : Let the system (1) be controlled by the M- 
MRAC scheme given by (5), ( 6 ), and (7). Then all closed-loop 
signals are ultimately bounded. 

Proof: Consider the following candidate Lyapunov func- 
tion 

V(t) = e T (t)Pe(t) — tr([0(f)0 T (f) + d(t)d T (t)] A) , (16) 

where d(t) = d.(t) — d. Computing its derivative along 
the trajectories of the systems (11) and (7), and using the 
properties of the projection operator, it is straightforward to 
obtain the inequality 

V(t) < —e T (t)(Q + 2XP)e(t) + 2e r (t)PB[d(t) — d] 

< -||e(f)|| [o||e(f)|| - 2d»] , (17) 

where we denote a = A m j n (Q) + 2A m i n (P)A and d * = 
\\PB m A(d(t) — rf) with A m i n (S l ) being the minimum 
eigenvalue of the matrix S. It follows that V(t) is negative 
semi-definite whenever a||e|| > 2d*, which along with the 
properties of the projection operator imply that the closed-loop 
error signals e(t),8(t), h(t) are uniformly ultimately bounded. 
The boundedness of e°(t) follows from the inequality (13). 
Since A rn is Hurwitz, the bounded signals r(t) and e(t) 
produce a bounded signal x m {f). Therefore, x(t) is bounded. 
It follows that u(t) and u(t) are bounded as well. ■ 


B. Transient behavior of the tracking error 

The projection operator in the adaptive laws (7) guarantees 
the inequalities 

\\m\\<e*, \\d(t)\\<d* , (18) 

therefore 

tr ([9(t)9 T (t) + d(t)d (i)]A) < (3 , 

where /3 = 49* 2 + 4 ft* 2 . From Theorem 5.1 it follows that 
V(t) < 0 whenever 

Arp Q 

V(t) >F* =A max (P)^ + ^. (19) 

cr 7 

Therefore, it follows that the trajectories stay inside the 
Lyapunov level set 

L={(e,M): V(e,8,d) = V*} . (20) 

From the definition of V (t) we have 

A min (P)||e(f )|| 2 < e T (t)Pe(t) < V(t) < C* . ( 21 ) 

Hence, the following conservative bound can be derived 


|e(f)|| £oo < c = 


c(P) 4 d 2 


P 


i(P) CL 2 jA min(P) ' 


( 22 ) 


Since the inequality (22) holds uniformly in t, the bound 
||e(f)|| £oo < c follows. We notice that the second term in 


the square root can be arbitrarily decreased by increasing the 
adaptation rate 7 . The first term is independent of 7 , but can be 
arbitrarily decreased by increasing A. This is not the case for 
the conventional MRAC design. Indeed, when A = 0, the only 
design parameter that affects that term is A m ; n (Q)- However, 
increasing A m i n (<3) scales also P, hence increases d*. 

The bound on ||e°(f)||£ oo follows from the inequality (13) 
and has the form 


||e 0 (f)|| Coo <c(l + Afc m ) (23) 

We notice that the derived bound on ||e°(f)||£ oo cannot be 
arbitrarily decreased by increasing the design parameters A and 
7 . If we set A = co^/ 7 , where the proportionality coefficient 
Co will be defined shortly, the following asymptotic bound can 
be written 


lim || e 

7 — >-oo 


I Coo 5) 


I Amax(P) 


d 2 + 


c oP 


. ^min(P) * ' A mi „(P) 
which can be decreased by increasing A m ; n (P). 


(24) 


C. Transient behavior of the control signal 

To investigate the behavior of the control signal with respect 
to design parameters we recall that u(t) does not explicitly 
depend on design parameters A and 7. Instead, u(t) depends 
on 7 through the adaptive laws, and u(t) depends on A through 
the tracking error dynamics. Therefore, for the purpose of this 
subsection we assume that r(£) has bounded time derivatives. 
This assumption is only needed for the analysis purposes and 
is conditioned on the way the bound on the control signal is 
derived. Alternatively one could use the integral representa- 
tion of the parameter estimates and the error signal e(t) to 
derive an integral equation for the u(t) without assuming the 
differentiability of r(t). However, for more transparency we 
use the differential form of the equations. 

Differentiating u(t) and substituting the adaptive laws we 
obtain 

u{t) = -7 [p{t)l q - H{f)]Bl l Pe{t) - r a (t) , (25) 

where we denote p(t) = (f) T + 1 , r a (t) = 9 T (t)(f)(t), 
and H(f) = G(d)(t> T (t)cj>(t) + G(d). We notice that the 
terms p(t), r a (t) and H(t) do not explicitly depend on the 
design parameters A and 7. Moreover, from the results of 
the previous subsection it follows that all signals involved 
in the equation (25) are bounded. In particular, there exist 
positive constants a\, < 22 , 0:3 such that ||p(£)||c cx) < a 1 , 
||p(£)||c 00 < «2 and ||r 0 (£)||/; oo < 03 . We also notice that the 
matrix F(t) = p(t)I p — H(t ) is symmetric and positive semi- 
definite, since ||G(0)|| < 1 , which follows from the properties 
of the projection operator. 

Differentiating the equation (25) with respect to time and 
using the equation (15) we obtain the following second order 
differential equation 

u(t) + Xu(t ) + 7 F(t)Lu(t) = -7 F{t)B] n PA m e{f) 

-7 F(t)Bj 2 Pe(t) + 7 F{t)Lu°(t) - A r a (t) - r a (t) ,(26) 



where L = B^P B m A. It follows from the definition of 
the projection operator that and Theorem 5.1 that F(t) is 
piecewise continuous and bounded. Since all the terms in 
equation (26) are bounded, it can be considered as a second 
order linear equation with time varying coefficients in u(t). 
Although the equation (26) is non-autonomous, it can be still 
inferred that the adaptation rate 7 determines the frequency 
of the control signal u(t). Therefore, increasing 7 increases 
the oscillations in u(t) as it is the case for the conventional 
MRAC design. On the other hand A determines the damping 
ratio. Therefore increasing A suppresses the oscillations in the 
control signal u(t). That is, by selecting a proper value for 
A the desired performance can be achieved. This is the main 
difference from the MRAC design, which results when A = 0. 

We select a proper value of A from the perspective of 
minimizing the norm bound on the control signal u(t) in 
transient. To this end we notice that selection of the initial 
parameter estimates inside the convex sets defined by the 
projection operator results in H(t) = 0 on some initial 
interval [0 t\\ by continuity. Therefore, F(t) = p(t)I p on 
[0 t\}. To simplify computations we notice that the matrix 
L is symmetric and positive definite, therefore there exists an 
orthogonal matrix T such that D = TLT t is diagonal with 
positive entries da, i = 1, . . . ,p. That is, introducing a new 
variable v = Tu (v° = Tu°), we can write the equation (26) 
in the following form 


v(t) + A v(t) + 7 Dp(t)v(t) = 

7 p(t)Dv°(t.) + 7 Zi(f) - A z 2 (t) - z 2 {t ) . (27) 

where 27(f) = T[-p(t)B^PA m - p(t)B^P]e(t) and 
22(f) = Tr a (t). Let po = 1 ^ 2 Ql . Then for each component of 
vector v(t) we can write the following equation 

Vi{t) + \vi{t) + 'ydup 0 Vi(t) = 7 p(t)duv°(t) 

+iz u {t) - Xz 2i (t) - i 2 (f) + 7[cfoPo - dup(t)]vi(t) .(28) 


Since p(t) > 1, this equation is in the form of equation (38) 
from Appendix with k(t) = dup(t), 2a = A, ojf = 'ydupo and 
three external inputs 


7 p(t.)duv°(t), 




Z2i(t) ■ 


Applying the inequality (48) we obtain 

|w(f)| < c 3 yj vf (0) + vf {Q)e~ vt + a\ ||u° (f ) || £oo (29) 

(gi + 2) v / po | 


+ -l-\\zu(t)\\c a 

da, 


/did 




for each i = 1 , ,p, if A > 2 yfUTi = y/jdupo. Therefore 
selecting 


A = sj2 1 d?(a 1 + l) 1 (30) 

where d° = max{cfii, i = 1 , . . . ,p} results in the inequality 

M0I < c 3 ^(0) + Vi(0)e~ vt + ai |k 9 (f) || a*, (31) 

1 ,, (Cl + 2 )y/~P0 , 


\\zii(t)\U 


+ 


\fdol 


\c 0 


for all components simultaneously, where do = ininfci,,. i = 
1, . . . ,p}. A similar inequality holds (as it is derived in the 
appendix) for the vector v(t) 

IKf)|| < C4e-" t + a 1 ||rt°(f)|| £oo (32) 

+ Uz 1 mc^ + {ci+ / p^ ro \\z 2 mc ao 

do V «o7 

Since T is orthogonal with ||T||2 = 1, we have ||u(f)|| = 
\\v(t)\\, ||z 2 (f)|| = ||r- a (f)|| and 

\\zi(t)\\ Coo < [arWBlPA^+^WBlPW] ||e(f)|| £oo . 


Combining the above relationships we obtain 


|«(t)|| < C4e-" t + Ql ||tx 0 (f)|U oo +c 5 ||e(f)|| £o 

(ci + “2)yfpo .. 

+ 7= lko(f)|Uoo , 

v«o7 


(33) 


where c 5 = ^ [a^B^P A m \\ + a 2 \\B^P \\] . 

It is easy to see that the bound on ||u(f)|| has the form 

C6 


u{t )\ | £oo < c 4 e v + a\\\u (f)||£ + — , (34) 


where 


I A n 


C6 = Co\ 




KJP) ' W py 


and the constants 01 . c\, cr, do not depend on the adaptation 
rate 7. We notice that the first term on the right hand side of 
(34) represents the exponentially decaying effect of the initial 
conditions w(0) and w(0), and the last term can be decreased 
by increasing 7. Therefore we can conclude that for large 7, 
the adaptive control u(t ) behaves similar to the ideal control 
u°(t). That is unlike the conventional MRAC design, the fast 
adaptation in the M-MRAC design does not generate high gain 
effect. 

The transient behavior of the error signals is summarized 
as follows. 

Theorem 5.2: Let the system (1) be controlled by the M- 
MRAC scheme given by (5), (6) and (7), and the design 
parameters are chosen according to equation A = co^/7 with 
co = i/2cf°(ai + 1). Then, the inequalities (22), (23) and (34) 
are true. 


VI. Simulation Results 

For the simulation we consider a wing rock motion of a 
slender delta wing given by the equation [14] 

4>{t) = /(</>, </>) + bu(t) + d(t) . 
where <j>(t) is the roll angle, 

/(</>, (j>) = ai4> + a 2 </> + ao\4> 2 \(j> + a 4 |</>|</> + a 3 </> 3 . (35) 

with ai = —0.0186, a 2 = 0.0152, a 3 = —0.0625, <24 = 
0.0095, <25 = 0.0215, 6 = 1. The disturbance d(t) rep- 
resents unknown atmospheric effects and is a square wave 
of amplitude 0.15 and of frequency 0.5 rad I sec. Only the 
sign of the control effectiveness b is assumed to be known 
(positive). The parameters of the reference model are chosen 



as follows A rn = [0 1; —1 — 1.6], b m = [0; 1]. The external 
input to follow is a step command of magnitude —15 degrees 
at t = 15. We run two simulations respectively from the 
small initial conditions (6 deg., 3deg/sec) and large initial 
conditions (30 deg., 10 deg/ sec). The adaptation rate is set 
to 7 = 10000 with A defined according to equation (30). 
The simulation results are displayed on Figures 1 and 2 
respectively. It can be seen that a good tracking is achieved 
for both output and control signals, and the later does not 
exhibit any high frequency oscillations even for the selected 
high adaptation rate. The disturbance effect is completely 
attenuated. 


Bank angle response in degrees 



Aileron deflection command in degrees 



Fig. 1. Tracking a step command from the small initial conditions. 


Bank angle response in degrees 



Fig. 2. Tracking a step command from the large initial conditions. 


VII. Conclusions 

We have presented design and performance analysis of M- 
MRAC architecture for a class of uncertain systems subject 


to bounded disturbances. It is shown that the systems’ input 
and output tracking errors can be decreased as desired by 
increasing the adaptation rate, when the error feedback gain 
is selected according to derived rule. This design method 
prevents high frequency oscillations in the control signal, 
which are unavoidable in conventional MRAC systems. The 
performance of M-MRAC is demonstrated on a benchmark 
problem of controlling wing rock motion of slender delta 
wings in a turbulent atmosphere. 

Appendix A 

An upper bound for a second order LTV SYSTEM 
Consider a second order time variant linear system 

x(t) + 2ax(t) + jk(t)x(t) = b\f(t) + b 2 f{t) (36) 

with x(0) = Xq, T(0) = Xq, where 7 > 0 is a constant 
parameter, k(t ) is continuous with k* > k(t) > fc* > 0 and 
has a bounded derivative. The function f(t) is assumed to be 
piecewise continuous and bounded. The equation (36) can be 
written in the matrix form as 


z(t) = Az(t ) +Bf(t ) 


(37) 


where 

z(t) = 


x(t) 

x(t) 


A = 


0 1 
—7 k(t) — 2 a 



We are interested in minimizing the upper bound on x{t) 
by the choice of the parameter a. To this end we introduce 
notations ui 2 = 'yko, ko = - ~^ k * , a = (ui and represent the 
system (37) in the following equivalent form 


z(t) = Dz(t ) + Bf(t) + C[w 2 — 'yk(t)]x(t) (38) 


where 


For the convenience of derivations we decompose z(t) into 
initial response z l (t) of the homogeneous system 

z*(t) = Az\t) (39) 


D = 


— w 


0 1 

2 -2 Cw 


with the initial condition zq = [cco io] T > an d force response 
z? (t) of the system (38) with zero initial conditions, which 
can be represented in the equivalent integral form 


= fo G ( t ~ r)Bf{r)dT + 
fg G(t — t)C[uj 2 — 7 k(t)]xf (r)dT . (40) 


Here, G(t ) = e Dt is the state transition matrix, which can be 
computed by direct integration (see for example [5]). 

In order to minimize the bound on z{ (t) = xJ (t) we 
compute the L\ norm of the elements in the first row of matrix 
G(t). For 1712(f) we obtain 




\\9l2(t)\\ Cl = 


+1 


2^3 — 1 


C<1 

C>1 


(41) 


where 8 = wy/|l — C 2 | and 0 < <p = tan 1 (^j) < f- 
Obviously, ||< 7 i 2 (t)||,c 1 reaches its minimum of A for all 
( > 1. On the other hand, the C\ norm of 511 (f) is computed 
to be 


IlimWIk = 



C<i 

C>i 


(42) 


It can be shown that the minimum of || 5 n(t )||£ 1 is reached 
at some £* < 1. The analytical computations are somewhat 
involved, but numerical computations results in CJ = 0.66 
with the minimum value of — , where c = 0.8026. 

Since there is no common minimum point for the entries 
of G(t), one can use a ’’suboptimal” value ( = 1, which is 
good enough for our purposes. In this case ||5n(i)|| ; c 1 = (j 
and || 512 (f) || £! = A- We notice that selecting a larger value 
of C while leaving ||5i 2 (f)||£ 1 intact, increases ||5n(f)||,c 1 
proportional to i. Therefore, we can select any > 1 with 
||5ii(f)IUi = — , where <7 > 2 is determined by the selected 
( and is independent of 10 . 

Next we compute Coo bound on x f (t) 


x f (t) = f* [61511 (f - t) + 62512 (f - r)Jf(r)dT 

+ So 512 (f - r)[w 2 - 7 k{t)]x f (T)d,T . (43) 

Since ||u ; 2 — / yk(t)\\c 00 = w 2 — 7 /c* = 7 fc ~ fc * , we obtain (see 
[7], p. 199 for details) 


Ik'WIUoo < ll/WIUJNIIffiiWIk + |ft2|||fl 1 2(f)lk 1 ] 

+ 7 ^||x/(f)||£ oo || 5 i 2 (f)||£ 1 . (44) 


Substituting the C\ norm values and solving the resulting 
inequality for ||a7 (f)|U 00 we obtain 


1 - 


k* - k t 


2ko 

which results in 


II^WIUcc < 


ci|6i| |6 2 | 


UJ 


UJ * 


ll/Wlkoo , (45) 


Uoo < 


ci 1 61 1 y/kp | 6 2 | 


\m\u 


(46) 


To obtain a bound for z*(f), we recall that according to 
Theorem 8.7 [13] the origin of the system (39) is uniformly 
exponentially stable, since ||A(f)|| is bounded, ||.A(f)|| is 
essentially bounded, and the point wise eigenvalues of matrix 
A(t) have negative right hand sides. Therefore 

||z < (f)||< C 3 ||*(0)||e^ (47) 


for some positive constants C 3 and v. 

Sice (49) is true for each component of z l (t), adding the 
corresponding inequalities we finally arrive at 


|*(f)|< C 3 ||z( 0 )||e-^ + 


ci|6i|v^ 

k*sjl 


1^1 
fc* 7. 


Il/Wlkoo (48) 


for all a > VT^o- 

We notice that when x and / are 5 -dimensional vectors 
in the equation (36), then z = [27 x\...x q x q ] T and the 
matrices A, B, C, D and G have repeated block structures. 


Therefore the equation (44) and the upper bound (46) hold 
for each component x{(t) of vector x? (t) with /(f) replaced 
with fi(t). On the other hand, the inequality (47) is true for 
the 25 -vector z(f), hence it is true for the vectors x^(t) and 
x/ (t). That is 

||* l (f)|| < c 3 e _l/ V Noll 2 + ||*o || 2 = c 4 e _l/t (49) 


It follows that the inequality 

\\x{t)\\< Ci e- vt 


ci\bi\\fkp j&2 
fc*7 J 


\m\\< 


(50) 


holds in the vector case as well, when a > 1 / 7 fco- 
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